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Abstract 

We review the Coulomb gauge in QCD and report some recent results. The 
minimal Coulomb gauge is defined and the fundamental modular region, a region 
without Gribov copies, is described. The Coulomb gauge action is expressed in 
terms of the dynamical degrees of freedom with an instantaneous Coulomb interac- 
tion, and its physical meaning is discussed. The local Coulomb gauge action with 
phase-space and ghost variables is derived, and its BRS-invariance and renormaliz- 
ability are reviewed. It is shown that the instantanteous part V(R) of g 2 Doo(R,t), 
the time-time component of the gluon propagator, is a renormalization-group in- 
variant V(R) = /(RAqcd)/R, and that the contribution of V(R) to the Wilson 
loop exponentiates. It is conjectured that V(R) ~ k cou \R at large R, and that 
K coui provides an order parameter for confinement of color even in the presence of 
dynamical quarks. 
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1 Introduction 



Perturbatively, in each order, the Coulomb gauge must agree with Lorentz-covariant 
gauges. So why bother with the Coulomb gauge? 

Use of the Coulomb gauge is motivated by the conjecture JIJ that non-perturbatively 
there is an instantaneous confining color-Coulomb potential V(r) which rises linearly 
V(r) ~ K c r at large r. Here the color-Coulomb potential is defined by 

gf 0) D$(x) = gf 0) (TA i o\x)A i o ) (0)) = V(r)8(t) + non-instantaneous, (1) 

where x M = (t,x), r = \x\, A^ = (Aq,A), and the index (0) refers to unrenormahzed 
quantities. In the Coulomb gauge A is transverse, V • A = 0. The quantity k c is a string 
tension, characteristic of the Coulomb gauge, that is finite in physical units. 

Most of the results which will be reported here are described in more detail in 0. 
Among other things, it is shown that V(R) is a renormalization-group invariant, and 
thus depends only on physical masses, V(r) = r~ l g1(rhQCD)- Here g c (rh.QCD) is a 
running coupling constant characteristic of the Coulomb gauge, which is independent of 
the renormalization mass \x and cutoff A, but which satisfies the usual renormalization- 
group equation. 



2 Basics and Notation 

The Yang-Mills field tensor is written 

^, = ^-94 + ^x4 (2) 

where (A^ x A u ) a = f abc A b ^A c v , and f abc are the structure constants of the Lie algebra of 
SU(N). The classical field equations are given by 

D^(A)F flu = j qui/} (3) 

where D^A) = 8^ + g^A^x is the gauge covariant derivative in the adjoint representa- 
tion, and the metric tensor has signature (— , +, +, +). Under a local gauge transforma- 
tion g{x) G SU(N), the connection A^ = A"M^ transforms according to 

9A^ = gA^ + gd^. (4) 

Its infinitesimal form, for g(x) = expo;(x) ~ 1 + u)(x), is given by 9 A^ = A^ — D^uj. 
The color- Gauss law 

D-E = p qu (5) 
is a constraint that will be essential to the confinement picture in the Coulomb gauge. 
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3 Minimal Coulomb gauge 



Gribov was the first to realize that the transversality condition does not uniquely fix the 
gauge in non-Abelian gauge theory 0. Thus, there exist different transverse configu- 
rations A and Ai, with V • A = V • A\ = 0, that are related by gauge transformation, 
A\ = 9 A. Here A and A\ are called "Gribov copies" of each other. As Gribov recognized, 
this has a significant effect on the dynamics because two configurations that look very 
different may in fact be Gribov copies of each other and thus represent the same physical 
configuration. The physical configuration space is the quotient space V = A/Q, which is 
the space of configurations A modulo local gauge transformations g. It is gauge-invariant 
and is also known to be topologically non-trivial (Singer's theorem) ||. A complete gauge 
fixing provides a concrete description of the quotient space by means of representative 
configurations. We now provide such a description by means of the minimal Coulomb 
gauge. 

We shall fix a gauge by minimizing a conveniently chosen functional on each "gauge 
orbit" [§, H, H, a gauge orbit being the set of all gauge transforms of one particular 
configuration A. ("All" must be suitably defined in the continuum theory. However the 
lattice analogs of quantities introduced below are well defined, as is the minimal lattice 
Coulomb gauge [0, §].) The space of gauge orbits is the quotient space A/Q. 

A "minimizing functional" for each time t is introduced, which is the restriction of 
the Hilbert norm to a generic gauge orbit, 

F x {g) = rA\\\ (6) 

where 

\\A\\ 2 = J d 3 x \A(t,x)\ 2 . (7) 

is the Hilbert norm of A at time t. To fix the gauge, we choose as representative for 
each gauge orbit that configuration for which F^(g) achieves its absolute minimum. The 
fundamental modular region (FMR) is the set A of representatives of all the gauge orbits 

A = {A: F x (l) < Fjfa) for aRg} . (8) 

This is a region without Gribov copies. Well, almost without Gribov copies. There 
are cases of "accidental" degeneracy, where the absolute minimum is not unique. The 
degenerate absolute minima are Gribov copies that our gauge-fixing procedure has not 
eliminated. One may show that they occur only on the boundary of A. These Gribov 
copies must be identified topologically for A to represent the quotient space. 



2 



This specifies a complete gauge fixing, but it is not very explicit. However from the 
existence of the minimizing functional, one may deduce explicit properties of the FMR. 
Any minimum, which may be relative or absolute, is a stationary point of F^(g), so 
5Fj(l)/8g = for all A e A. A simple calculation gives 5Fj(l)/5g = -2V • A, and 
we conclude that the Coulomb gauge condition V ■ A = is satisfied for all A e A. 
In addition there is a second condition which holds namely that the matrix of second 
derivatives of the minimizing function is non-negative at an absolute or relative minimum, 

— * — * — * — * 

and a simple calculation gives — D ■ V > 0, where D = D{A). Here the Faddeev-Popov 
operator M(A) = —D ■ V has made its appearance. It is symmetric for A transverse, 
—D • V = —V • D. These two conditions define the Gribov region Q. It is the set of 
relative and absolute minima of the minimizing function, 

n = {A : V • A = and - D ■ V > 0}, (9) 

that contains the fundamental modular region A as a proper subset A C fl. One may 
show that A and Q are both convex, that both are bounded in every direction (this is 
proven below), and that A and fl have some boundary points in common. 

4 Color- Gauss law 

We have remarked that Gauss's law D ■ E = p qu is essential for the confinement picture 
in the Coulomb gauge. To solve this constraint we pose E = E tT — V$, where $ is the 

— * — > 

color-Coulomb potential, and V • E tr = 0. Gauss's law reads 

-D-V$ = p. (10) 

Here we have introduced a color charge density 

P = -9(0) A x Ef + p qu (11) 

which includes only that part of the gluon charge density that comes from E tT . Gauss's 
law fixes $, 

$=(-D-VrV (12) 

This equation, which appears neither gauge nor Lorentz covariant, in fact holds in every 
gauge and every Lorentz frame. However in the minimal Coulomb gauge, A is transverse, 
A = A tv , so p depends only on the dynamical degrees of freedom of the gluon A tr and 
E tT . 
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Moreover, as we have seen, the Faddeev-Popov operator in the minimal Coulomb 
gauge, M(A) = —D ■ V is strictly positive in the interior of the fundamental modular 
region, M(A) > 0, so the inverse which appears in $ = M~ 1 (A)p is well defined, except 
possibly for some points on the boundary of A where M(A) may have a zero eigenvalue]] 
Thus in the minimal Coulomb gauge there is a coherence between the gauge fixing and 
the dynamics expressed by the Gauss law constraint, in the sense that this constraint 
may be solved by inverting a positive operator. 

The singularity of M~ 1 (A) strongly enhances the color-Coulomb potential $ as parts 
of the boundary of A are approached. This is an expression of the famous anti-screening 
property of QCD and, as we conjecture, of confinement. Indeed, the Green function 
which solves the color Gauss law constraint is given by 

M-\A) = (-V 2 - <7 (0) A x V,)- 1 (13) 

The first term —V 2 is strictly positive. However the second term has zero trace on the 
color indices, and thus has positive and negative eigenvalues. For some subspace there is 
a cancellation between the first and second term, which gives a small denominator. 

To see this in detail, consider the ray in A-space defined by XA, where A > is a 
positive parameter, and A ^ is arbitrary but fixed [[|. For sufficiently small A, the 
Faddeev-Popov operator M(XA) = — V 2 — Ag(o)^4j x V« is strictly positive. As A increases, 
the lowest (non-trivial) eigenvalue of M(XA) decreases toward 0, giving an enhanced 
color-Coulomb potential, as we now show. For let wbea vector for which the expectation 
value of the second term of M(A) is negative, E = (u, — <7(o)Aj x V; w) < 0. Such a 
vector certainly exists because this operator has zero trace. We have (u, M(XA) u) = 
||Vo;|| 2 + XE, with E < 0. This expression starts out positive and decreases linearly 
as A increases, eventually becoming negative. So for A sufficiently large and positive, 
the lowest (non-trivial) eigenvalue hq(XA) of M(XA) is certainly negative. Consequently 
fio(XA) is a continuous function of A that starts out positive and goes negative, and thus 
it inevitably passes through 0. (One may show that the decrease is monotonic.) When 
Ho(XA) = 0, the point XA lies on the boundary dQ of the Gribov region Q. The color- 
Coulomb potential $ = M _1 p is strongly enhanced as /j,q(XA) — > 0. Incidentally, the 
above argument also shows that the FMR is bounded in every direction A. 



1 Actually, M(A) has a trivial null space consisting of constant functions, Vw = 0, and the inverse 
M _1 is taken on the orthogonal subspace. Only V<& appears in E = E tr — V$, so the arbitrariness of 
$ under $ — * $ + oj is of no consequence. 



4 



5 Coulomb-gauge action in terms of dynamical de- 
grees of freedom 



A serious obstacle to the use of the Coulomb gauge in QCD has been doubt about 
whether the theory is renormalizable and thus whether the gauge really exists at all. In 
this section we indicate the nature of the problem and explain the solution physically by 
integrating out the constraints so only the dynamical degrees of freedom remain. 

The Faddeev-Popov formula for the Euclidean functional integral in the Coulomb 
gauge reads, 

Z = ( d A A5iy-A) det(-V--D) exp(-S d ), (14) 

J A 

where 

S d = (l/2) Jd'x(E 2 + B 2 ), (15) 

Ef = -Fqj, and Bf = Fj k for k cyclic. The integration over A is restricted to the 
fundamental modular region, A, but we temporarily ignore this point and discuss the 
perturbative expansion which is insensitive to this restriction. 

The above expression certainly cannot be perturbatively renormalized in the usual 
way, as one may see by considering the closed ghost loops that represent the expansion 
of the Faddeev-Popov determinent. The free propagator of a ghost with momentum k is 
1/k 2 . It is independent of k and corresponds to the instantaneous space-time propagator 
S(t)/r, so a closed ghost loop is an instantaneous closed loop in the Coulomb gauge. The 
insertion of the simplest ghost loop into a gluon propagator of momentum p gives the 
loop integral 

W^svbp < 16 > 

The ko integration diverges because the integrand is independent of ko- This occurs in any 
number of spatial dimensions, the so the divergence cannot be regularized by dimensional 
regularization. Moreover the divergence is non-polynomial in p. One suspects that this 
divergence is cancelled by other diagrams. The problem is to make this cancellation 
manifest . 

To do this we introduce the Gaussian identity 

1 = J d 3 Uexp[- J d 4 x {l/2){U + iE) 2 ] (17) 

into Z, which gives 

Z = [ dA tr d 3 TldA det(-V--D) exp(-S'i), (18) 

J A 
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where 

Si = f d*x [iU - E + (1/2) (n 2 + B 2 )}. (19) 

Here we have used the constraint expressed by 5(V - A), to integrate out the longitudinal 
part of A, so A = A tT in the integrand, and the integation over A tT remains. 

With II • E = II • (dA tr /dt — DAq), we observe that the action is linear in Aq. 
Integration over Aq imposes the Gauss law constraint, 5(D-U). We separate the tranverse 
and longitudinal parts of n according to IT = n tr — V^, where $ is the color-Coulomb 
potential, and we have and d 3 Tl = constdII tr The Gauss law constraint reads 5(M<& — 
p), where M = — D ■ V is the Faddeev- Popov operator, and 

p = -g (0) Af x nf . (20) 

We next effect the integration on $ using 

J d<S>detM 5(M$ - p) /($) = f(M~ l p). (21) 

This absorbs the Faddeev- Popov determinent, and thereby achieves the goal of eliminat- 
ing the unwanted non-renormalizable instantaneous closed ghost loops. 

We obtain finally the partition function expressed as a functional integral over the 
dynamical phase-space variables, 

Z = [ dA tT dU tT exp(-S couI ), (22) 

where the Coulomb-gauge action is given by 

,S coul = J d 4 x {ffl tr • dA tr /dt + (l/2)[(n tr ) 2 + (VM"V) 2 + B 2 )]}. (23) 

This formula is physically transparent. The dynamical variables are the canonically 
conjugate variables A tT and n tr that remain after the constraints are integrated out. 
This formula could also have been obtained by the method of canonical quantization of 
constrained systems. As desired, there are no instantaneous closed loops because the 
Faddeev-Popov determinent is gone. There remains the instantaneous color-Coulomb 
interaction 

5 inst = J d 4 x (VM^V) 2 = / dtd 3 yd 3 z p(t, y) Q(y, z; l tr (t)) p(t, z), (24) 
with kernel 

Q = M^i-V^M- 1 . (25) 
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This kernel contains the inverse Faddeev-Popov operator M _1 (twice). Recall that 
the lowest eigenvalue of M approaches as a generic point on the boundary dQ of the 
Gribov region Q is approached. Consequently the Coulomb-gauge action S cou i diverges 
as dQ is approached, and exp(— Seoul), regarded as a function of g^A tr , is strongly 
suppressed by an essential singularity on dQ. The perturbative expansion in powers of 
#(o)^4 tr reflects this strong suppression. Therefore, although the perturbative expansion 
is insensitive to the restriction to the fundamental modular region A C £1, one may 
speculate whether the dynamical suppression by exp(— S cou i) as dQ is approached, is 
sufficient for the perturbative expansion to capture in some sense the non-perturbative 
features of the theory, for example by resummation. 



6 Renormalizability of the local Coulomb- gauge ac- 
tion 

Formula (|22| ) is physically transparent, but it is unsuitable for renormalization which 
requires a local action. We therefore start over and express the Faddeev-Popov formula 
(|14D in terms of the local Faddeev-Popov action 

Z = J d^AdCdCdXexpi-S-Fp), (26) 

where 

S FP = J d 4 x{ (1/2) {E 2 + B 2 ) + iVC ■ DC + iVX ■ A]. (27) 

Here C and C are the usual Faddeev-Popov ghosts and A is a Lagrange multiplier that 
enforces the gauge constraint. The instantaneous closed fermi-ghost loops have been 
reintroduced, but we have seen that cancellation of closed instantaneous loops is simpler 
in the phase space representation. We therefore insert the Gaussian identity eq. ( |TTD 
into the last expression and obtain 

Z = J d*Ad 3 IldCdCd\exp{-S), (28) 

where 

S = J d A x[ iU-E+ (l/2)(n 2 + B 2 ) + iVC -DC + tVX-A }. (29) 

We shall show that the instantaneous closed fermi ghost loops cancel against the instan- 
taneous closed bose loops. 
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Before doing so, we discuss some properties of this action. The equations of motion 
of A and II read 

D-U=0 

U = E. (30) 

Thus in the local formulation of the Coulomb gauge, Gauss's law is satisfied as an equation 
of motion, despite the coupling to ghosts. This is the underlying reason why g^A^ = 
gr r \Ay is invariant under renormalization. 

Consider the BRS transformation defined by 

sA, = -DpC, sC = -(1/% (0) C x C 
sC = X s\ = 

sUi = C x IT;. (31) 

Here A, C, C, A transform in the usual way, and II transforms covariantly, so s is nilpotent, 
s 2 = 0. Because s effects an infinitesimal gauge transformation on A^, the Yang-Mills 
fields also transform covariantly, sEi = C x Ei and sBi = C x B{. The local action S 
may be written 

S = J d A x[ (1/2)(II + xEf + (1/2)(E 2 + B 2 ) + s{iVC -A)}, (32) 

which makes it obvious that S is BRS invariant, sS — 0. 

The only place in which a time derivative appears on the local action (|29|) is in the 
term ill ■ E = ill ■ (dA/dt — DA ), so from Noether's theorem, the BRS charge density 
is Pbrs = — • DC, and the conserved BRS charge is simply 

Qbrs = i / d 3 xCD • II = % J d 3 xC5S/5A , (33) 

which is intimately connected to Gauss's law. Conservation of this charge leads to a 
Ward identity which is specific to the Coulomb gauge, 

^(77-777- + -) = d o f ^(Ctt- + -)• (34) 



^ApSKp ' "J v 6A 

Integration of this equation over all time annihilates the right hand side and gives the 
familiar Zinn- Justin equation. So this identity is more detailed than the Zinn- Justin 
equation. It expresses the constraint imposed on renormalization by the color-Gauss 
law. From it one may deduce 

g m Af ] = g (r) At\ (35) 
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or in other words that goAo is a renormalization-group invariant 0. It follows that the 
time-time component of the gluon propagator g? -,D^(R,t) is also a renormalization- 
group invariant, as is its instantaneous part V(R). Thus it depends only on physical 
masses, V(R) = f(A Q c D R)/R- 

The breaking of Lorentz invariance by the Coulomb gauge-fixing occurs only in the 
BRS-exact piece of the action s(iVC ■ A). This allows one to control the breaking of 
Lorentz invariance, and one may show that the expectation value of a gauge-invariant 
quantity is Lorentz invariant. 

There is no difficulty in deriving Feynman rules from the local phase-space action 
S. For example the H — IT,- propagator is (k 2 5ij — kikj)/k 2 , the Hi — Aj propagator is 
ko(Sij — kikj/k 2 )/k 2 , and the remaining propagators are as usual. 

We now show that the the instantaneous closed fermi loops cancel exactly against 
instantanteous closed bose loops in every order of perturbation theory. The instantaneous 
closed fermi loops come from the term in the action iVC ■ DC. Now consider the term 
in the action 

m- E = 2(n tr -V$) ■ (OA/ dt- DA ) } (36) 

where we have separated out the transverse and longitudinal parts of II. The term 
zV$ • DA has precisely the same form as the term as the ghost term iWC ■ DC. For 
this reason it produces instantaneous closed bose loops in every order of perturbation 
theory that are equal, apart from sign, to the instantaneous closed fermi loops. Thus the 
non-renormalizable instantaneous closed loops cancel pairwise. Indeed this is why the 
ghosts are necessary. 

The last argument is not rigorous because it involves a cancellation of infinite quan- 
tities. However it may be made rigorous by use of the interpolating gauge, defined by 
the gauge condition — adoA + V • A = 0. Here a > is a gauge parameter, and in the 
limit a — > 0, the interpolating gauge approaches the Coulomb gauge. For finite a, the 
theory is renormalizable. One may show |I| that the cancellation of the unwanted fermi 
and bose loops occurs at finite a, modulo terms that vanish with a. This establishes 
renormalizability of the Coulomb gauge. 
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7 Exponentiation of the static color- Coulomb poten- 
tial 



Consider a rectangular Wilson loop 

W = (tr P exp( J g {0) A^W)) (37) 

oriented along the time axis, of dimension T x R. The anti-hermitian matrices A a repre- 
sent the SU(N) generators in some representation D. In the Coulomb gauge, the time- 
time component of the gluon propagator g^D^ has an instantaneous piece V(r)S(t). 
We shall evaluate the contribution W inst to the Wilson loop that comes from keeping 
only the instantaneous part of the gluon propagator, without attempting to evaluate the 
remainder. 

The instantaneous prt of the gluon propagator contributes all possible ladder diagrams 
to the Wilson loop. A typical ladder diagram, with time ordering, 0,ti < t 2 < ... < T, 
contributes 

tr(...A a2 A ai A ai A a2 ...) / dt 1 dt 2 ...V{R)V{R)... . (38) 

Jo<t!<t 2 <...<T 

The path ordering for a ladder diagram produces the nested product of A-matrices that 
is shown. The identity A a A a = —Co, where Co is the Casimir for the representation D, 
may be applied successively to the nested pairs of contracted A-matrices, 

tr(...A a2 A ai A ai A a2 ...) =tr(...)(-Cf,), (39) 

which allows us to evaluate all the traces. The contribution of instantaneous propagators 
to the Wilson loop exponentiates, giving 

W inst = (N 2 - 1) exp[-TC D V(R)/2]. (40) 

8 Order parameter for color confinement 

Recall that V(R) depends only on the physical mass Aqcd- This result, and the ex- 
ponentiation formula for W^ nst , inspire the conjecture that Wi nst obeys an area law, or 
equivalently that V(R) rises linearly at large R, 

V(R) ~ k cou1 R. (41) 

Here k cou \ is a string tension that is characteristic of the Coulomb gauge, but finite in 
physical units. 
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It is believed, and confirmed by numerical simulation in lattice gauge theory, that in 
gluodynamics the Wilson loop in the fundamental representation obeys an area law W ~ 
exp(-nTR), but that a perimeter law holds in the adjoint representation. Consequently, 
if eq. (pll)) is true, it cannot possibly also be true that W ~ Winst in all representations 
D, for eq. (f|l]) gives an area law for Wi nst whenever Co is non-zero. Nor does the Wilson 
loop W obey an area law when dynamical quarks are present, whereas we conjecture 
that in this case also V(R) rises linearly at large R. This is as it should be, because the 
color-Coulomb potential V(R) should not give the energy of a state which contains a pair 
of external quarks, for that depends on whether mesons are formed. Indeed the linear 
rise of V(R) at large R makes it energetically favorable in the presence of an external 
quark for a dynamical anti-quark to be produced from the vacuum, to form a meson, in 
which case the energy of a state containing an external quark pair at large separation 
does not rise linearly with the separation. Rather the color-Coulomb potential V(R) 
couples universally to color-charge. Thus if V(R) does rise linearly at large R, then 

k cou1 = lim V(R)/R (42) 

serves as an order parameter for color-confinement, which is otherwise lacking. 

It appears that this conjecture is susceptible of investigation by numerical simulation 
in lattice gauge theory. Configurations that are taken from the Wilson ensemble may be 
gauge fixed numerically to the minimal Coulomb gauge, and V(R) may then be found 
as the expectation- value of instantaneous (t = 0) part of g 2 D 0Q (R,t). 
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